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Interaction induced trapping and pulsed emission of a magneti- 
. . xally insensitive Bose-Einstein Condensate 
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PACS 03 . 75 . Mn - Multicomponent condensates; spinor condensates 
PACS 34 . 50 . Cx - Elastic; ultracold collisions 
PACS 67 . 85 . -d - Ultracold gases, trapped gases 

Abstract. - We demonstrate that atoms in magnetically insensitive hyperfine states (m=0) can 
be trapped efficiently by a Bose-Einstein Condensate of the same atomic species occupying a 
different hyperfine state. The latter is trapped magnetically. Hyperfine state changing collisions, 
and therefore loss of the trapped (m=0) atoms, are shown to be strongly inhibited in case of a 
low density of the confined atomic cloud. We monitor the transition from a 'soft' to a 'hard' 
effective potential by studying the backaction of the trapped (m=0) atoms onto the condensate 
which provides their confinement. The controlled outcoupling of the trapped atoms by shaping the 
condensate's wavefunction is explored. We observe a pulsed emission of atoms from the trapping 
region reminiscent of an atom laser. 
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Trapped ultracold atomic gases offer outstanding possi- 
bilities to model complex quantum systems. A paradigm 
is the Mott-Insulator phase transition demonstrated with 
atoms trapped in a lattice potential [1]. Thus, a key 
ingredient for advancing the possibilities to explore the 
quantum dynamics of many body ensembles with ultra- 
cold atoms is to improve existing methods or to find novel 
ways for the control of the external motion of atoms. Usu- 
ally static or time-dependent electro-magnetic fields are 
employed for this task [2-6]. A prominent example are 
optical lattices formed by counter propagating light waves 
who facilitated the observation of effects like the above- 
mentioned Mott-Insulator phase transition [1], second- 
order tunneling [7] or Joscphson oscillations [8]. In a 
similar fashion so-called radio- frequency dressed adiabatic 
potentials that emerge from a combination of static and 
oscillating magnetic fields [9, 10] have been used in order 
to coherently manipulate matterwaves [11] and to inves- 
tigate the decoherence dynamics of one-dimensional Bose 
gases [12]. 

Magnetic trapping relies on the coupling of the total 
angular momentum F to the magnetic field vector. The 
resulting potential is proportional to the magnetic projec- 



tion quantum number m ranging from — F to F. Hence, 
in the linear Zeeman regime the m = state is insensitive 
to the magnetic field and can therefore not be trapped 
magnetically. By contrast, in optical traps it is possi- 
ble to trap the m = component and to achieve atomic 
spinor Bosc-Einstcin Condensates (BEC) where the dif- 
ferent spin components can exchange population coher- 
ently [13]. This spin mixing dynamics has been studied 
theoretically [14,15] and experimentally, e.g. for 23 Na [16] 
and 87 Rb [17]. 

In this letter we show that the m = component can 
nevertheless be confined within a magnetic trap. However, 
the confinement is not provided by a potential due to the 
magnetic field but rather by a BEC formed by atoms in a 
state with m ^ 0. We show that for a large BEC density 
in the m/0 state the atoms in the m = state behave as 
if they were confined by a potential which is constituted 
by the condensate density. Eventually we demonstrate 
how atoms in the m = state are released in a controlled 
manner from the 'trap' if the BEC density is varied. This 
method complements current methods to outcouple atoms 
from a BEC, e.g. radio frequency [18] or Raman transi- 
tions [19] and facilitates the creation of an atom laser for 
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m = atoms. 

In the following we consider atoms in a F = 1 hyper- 
fine manifold and restrict ourselves to an effective one- 
dimensional description accounting only for the longitudi- 
nal dynamics. The transversal confinement may be pro- 
vided by an isotropic harmonic potential with a sufficiently 
high trap frequency u)±_ (associated with the harmonic os- 
cillator length ajj such that the transversal dynamics is 
frozen out (see chap. 1 in [20] and ref. [21,22]) i.e. only the 
corresponding ground state is occupied. The longitudinal 
motion takes place in ^-direction and the corresponding 
mean field equations for the evolution of the spinor com- 
ponents ^ m = ty m {x,t) are given by 
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with the atomic mass M, the total atomic density n = 
EL-i I*"/' G z = |*+i| 2 - |*-i| 2 and G+ = G*_ = 
V2 [^i^o + ^o^-i]- The coupling constants are go = 
2h 2 (a + 2a 2 )/{3Ma 2 L ) and Sl = 2h 2 (a 2 - a )/(3Ma 2 L ) 
where a / 2 are the s-wave scattering lengths for the 
scattering-channels with total spin and 2. The parti- 
cle number in the respective hyperfine state is calculated 
according to N m = J dx \ ^! m \ ■ ln this work we consider 
23 Na with the respective scattering lengths ao = 2.43 nm 
and a 2 = 2.75 nm. The axial trap frequency is chosen 
lu± = 2 kHz which gives rise to a transverse oscillator 
length a± = 2.97 x 10 3 nm. 

The spin-dependent potential under consideration is of 
the form V m (x) = mV(x) with V(x) forming a double- 
well. As shown in Refs. [9, 23] such a potential can be 
created using a standard magnetic trap of Ioffe-Pritchard 
type which is dressed by a homogeneous radio-frequency 
field. Instead of using the actual form of V(x) which is 
provided in Ref. [9] we rather model the double-well as 
V(x) = 7 + jycxp (-x 2 /a 2 ) + (l/2)Mu¥. This model 
potential, which is easily implemented numerically, covers 
the important features of the exact potential and can be 
made to match the experimental situation by tuning the 
parameters 7, 77, a and u>. For our calculations we choose 
7 = 0, 17 = 1.05 x 10 30 J, a = 28 Mm and uj = 71 Hz. A 
sketch of this potential is shown in fig. 1. While atoms 
in the m = 1 state are subjected to an overall confining 
double well potential, atoms in the state m = — 1 can be 
trapped only temporarily in metastable states of a single 
well which is unbounded from below. We assume this state 
to be stable over the timescale of interest. Additionally, 
the atom numbers and well separations are chosen such 
that there is negligible initial overlap between the three 
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Fig. 1: Preparation of the initial state: BECs prepared in the 
two wells of the m = 1 manifold and the single well of the 
m = — 1 potential. All m = — 1 atoms are transferred by a 
Raman transition to the m = state. 



atomic clouds located in the three wells. The initial state 
for our investigations is created by transferring the entire 
population of the m = — 1 into the m = state which 
experiences no confining potential. Such a transfer can be 
achieved experimentally by two sufficiently broad banded 
microwave or laser pulses providing a Raman transition 
via intermediate excited states. 

We prepare our initial state by relaxation, i.e. imagi- 
nary time propagation of a trial wave function. Since there 
is no initial overlap between the wavepackets each of them 
can be prepared independently We therefore introduce 
the number of particles confined in the left and right well 
by N lL = J^l^il 2 and N lR = f °°|*i| 2 , respectively 
(sec fig. 1). They obey N 1L + N 1R = N x . After hav- 
ing obtained the ground states in each potential well we 
set "Jo = ^-1 and subsequently put = assuming 
a perfect population transfer. Then the system is evolved 
in time according to eqs. (1) using the Adams-Bashforth- 
Moulton predictor-corrector method. 

Throughout this work we will consider the regime of low 
densities of the m = atoms and a small maximum over- 
lap region of the m = and m = 1 wavepackets in which 
case dtN-i becomes negligibly small. For analytical con- 
siderations we therefore make the very good approxima- 
tion vff-i = for all times (see below for a corresponding 
discussion) and arrive at the two equations of motion: 
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i.e. the m = — 1 component is not involved in the time- 
evolution at all. From eq. (2) we observe that, although 
there is no external trap acting on we can identify the 
term V e ff(x) = (go + gi)\^i\ 2 as an effective potential. 
However, this potential depends explicitly on the m = 1 
density and implicitly on the m = density as seen from 
eq. (3). In order to demonstrate that V e ff{x) can provide 
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Fig. 2: a: Density of the individual m components (m = 1 
dashed line, m = solid line) of the initial state with Ni = 10 
(Nm = Nil) and iVo = 10 3 . b: Density snapshot after the 
time t = 13 ms. c: Time evolution of the density of the m — 
component, d: Relative change of the number of particles 
A-ZV™(t) = JV "' ( ^~^"' (0) for JVo(0) = 10 3 . The change of the 

1 and m — — 1 components are 
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particle number of the m 



equal (upper lines). The lower line denotes the change of the 
particle number of the m — component. 



trapping we prepare an initial state with Nm = Nil = 
(1/2) x 10 5 and = 10 3 as shown in fig. 2a. During 
the first several milliseconds the non-stationary m = 
wave packet broadens until parts of it hit the atoms in the 
m = 1 state and \&i and Wo overlap. The situation at this 
instant of time is depicted in fig. 2b. The m = atoms are 
reflected back completely and are thus eventually trapped 
as displayed in fig. 2c, where the m = density over a time 
interval of 100 ms is shown. As can be seen the atoms are 
effectively confined to the interval — 50/Ltm < x < 50 /zm. 
Eq. (2) shows that the trapping is rooted in the density- 
density interaction between atoms in different spin states. 
A similar behavior is therefore expected for 87 Rb and for 
any other atomic species whose scattering lengths obey 
.9o + ,9i > 0. Moreover, the density-density interaction 
implies that the demonstrated trapping effect should also 
be observable if two different atomic species are employed 
instead of two different spin states of the same atom. This, 
of course, requires that the atomic density of one of the 
species is shaped accordingly. 

Let us now return to the question of particle loss due to 
the scattering process (m, m!) = (0,0) — ► (1,-1). In fig. 
2d we show the temporal evolution of the relative change of 
the particle number AN m (t) = j^&i!l^rai2J m each of the 
three m-channels for Nq(0) 



N (0) 

10 3 . The relative particle 
loss from the trapped atoms in the m = state is negli- 
gibly small for the time interval 100 ms and therefore eqs. 
(2,3) provide an accurate approximation. The scattering 
process (0, 0) — > (1, —1), and thus the loss of m = atoms, 
is strongly suppressed. For vanishing m = 1 and m = — 1 
components eqs. (1) reduce to the Gross-Pitaevskii equa- 
tion for the m = component only and spin exchange 
processes are therefore absent. The microscopic spin ex- 
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Fig. 3: Relative energy change AE(t) = 

m = component for No = 10 3 , Ni = 
N 1R = N 1L . 
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change process can therefore in the mean field picture only 
occur if there is an overlap of at least two different spin 
components. In our situation this overlap is negligibly 
small. In order to see this more quantitatively we derive 
an upper bound for the population rate of the m = — 1 
component on the basis of eq. (1). The particle num- 
ber obeys 6W_i = / dx [d t ^- X ^*_i + St*li*-i] • 
For short times the evolution of is governed by 

dt'fy-i = — i^-^o^i- Assuming the r.h.s. to be con- 
stant leads then to ^^i(t) w dt^-i t. Considering ex- 
clusively the process that populates $_iwc find dtN-i w 
/ dxdt® -idt®*^ w (gi/h) 2 f da;|*o| 4 |*i| 2 t. Introduc- 
ing a typical overlap region Ax between the m = and 
rn = 1 wave functions and the corresponding maximal val- 
ues for the densities Tio.max; ^l.max 

inside this region, we 
find the following upper bound for the particle number 
increase in the m = — 1 mode: 



6W_i < (-^-) 2 Axno max ni, max t. 



(4) 



Estimating the increase of the particle number in the m = 
— 1 according to eq. (4) we find with no max = 25 fitn , 
ni.max = 100 fim and an overlap of Ax ~ 20 fim a 
rate of dtN-i < 0.02t /(ms) 2 which is consistent with our 
numerical data. Beyond this estimate fig. 2d shows a step- 
like increase of iV_i which is due to the time-dependence 
of the overlap of an d ^i- 

Let us now address the question of a backaction of the 
ni = atoms onto the m = 1 atoms. As demonstrated 
above the effective potential V e ff(x) = (.9o + .9i) |*i| 2 can 
grant efficient confinement of the m = atoms. How- 
ever, according to the coupled eqs. (2,3) does not 
evolve independently of ^>o- Hence it is expected that 
the trapped atoms act back on the atoms in the m = 1 
state and thereby modify the trapping potential itself. 
In order to quantify this effect we calculate the effective 
energy E e ff - 

2a 2 h 2 
Map, 
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|*o| 2 |*i| 2 ) which is the energy contained in the 
component itself plus a contribution arising due 
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to the interaction of the to = with the m = 1 density. 

If l^i | 2 was static, i.e. it did not change its shape as 

a function of time, E e ff — E e ff{t) would be conserved. 

Conversely, energy exchange between the spin components 

will be reflected in a variation of E e ff over time. In 

fig. 3 we illustrate the time-dependence of the quantity 
A £( t ) = E cff m-Eef f (o) for No = 1q3 and twQ different 

particlc numbers iVi = 10 5 and 5 x 10 4 . The shape of the 
initial state is the same as discussed before, i.e. symmet- 
ric occupation of the to = 1 double-well potential (see fig. 
2a). We observe that for N\ = 10 5 the relative energy 
change is of the order of 5 % over the shown time interval 
of 10 ms. Here the effective potential can be approxi- 
mately considered as static. The situation changes if N\ 
is lowered. More concrete, for Ni = 5 x 10 4 we observe 
large oscillations of AE with peak values up to 35 % and 
the effective potential picture breaks down. 

So far we have been focusing on an initial state with a 
symmetric occupation of the m = 1 double well. We will 
now investigate the situation Nil ^ Nm where Nm is 
not large enough to provide a complete confinement of the 
m = atoms. In fig. 4a the evolution of the to = density 
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Fig. 4: a,b: The time evolution of the m — density for 150 
ms and an initial state with N = 10 3 , N 1L = 10 5 , N 1R = 
(a) and Nut = 7 x 10 3 (b). Atoms are released from the trap 
region and escape to the x > halfspace. c: Time evolution 
of the relative number of escaped m = atoms for (from top 
to bottom) N 1R = {0, 10 3 , 2 x 10 3 , 7 x 10 3 }. d: Number of 
trapped atoms N r — No — Nq 30 after t = 150 ms as function 
of N 1R for No = 10 4 (+), N = 5 x 10 4 (x) and N = 10 3 (A). 
The lines were generated by fitting iVo r = k\/Nir. 

is shown for the particle numbers No = 10 3 , Nil — 10 5 
and Nm = 0. In the case, where there is no occupation of 
the right hand well, the atoms are free to leave the trap- 
ping region to the x > halfspace. Escape to the opposite 
direction is prevented by the m = 1 wave packet in the left 
well (Nil) and thus atoms which are initially going to the 
left are reflected. They interfere with the atoms initially 
going to the right thereby giving rise to an emission of a 
sequence of distinct wave packets. In case of an initial oc- 
cupation of the right hand well, i.e. Nm ^ the situation 
becomes more complex. An example is shown in fig. 4b 



where Nm = 7 x 10 3 . Here Nm is not sufficiently large in 
order to provide confinement for all to = atoms, and a 
fraction of them escape to the x > halfspace. This emis- 
sion recedes drastically as soon as a sufficient number of 
atoms has escaped and the number of to = atoms in the 
trapping region has become so small that confinement can 
eventually be granted by the to = 1 atoms in the right well. 
This behavior is studied in more detail in fig. 4c where we 
plot the time evolution of the number of escaped atoms 
JVg sc = \^o\ 2 for different values of the atom num- 
ber Nm keeping Nq = 10 3 and Nil = 10 5 constant. In 
the case at hand wc define xq = 80 /im thereby ensuring 
that the considered atoms are definitely not trapped any 
longer since the density of the m = 1 component beyond 
Xo is sufficiently small. The number of escaped atoms 
increases monotonously as time passes. For t < 60 ms pe- 
riods of a steep increase are followed by plateaus where 
almost no emission is observed. During the latter time 
intervals those atoms which are reflected from the to = 1 
wave packet localized in the left well destructively inter- 
fere with the atoms that initially went to the right, caus- 
ing the emission to cease for a while. This is reflected in 
the pulsed release of to = atoms from the trapping re- 
gion that is clearly seen in figs. 4a and b. This situation 
is strongly reminiscent of a pulsed atom laser. For large 
times Nq sc /Nq saturates since the number of remaining 
atoms in the m = state is small enough to be confined. 
The number of confined atoms Nq t = No — Nq sc after 
t = 150 ms is plotted in fig 4d as a function of the occupa- 
tion number of the right well for different No. As expected 
the final number of trapped atoms increases with increas- 
ing Nir. The energy Eq of the to = component is ap- 
proximately proportional to the square of the occupation 
number No in case of the nonlinear term being dominant. 
Furthermore, if one regards the m = 1 component in the 
right well to act as a potential then the height Vm of this 
potential is proportional to Nir. The energy Eq oc Nq of 
trapped m = atoms should be less than the height of the 
potential Vm oc Nm created by the to = 1 atoms. This 
yields the scaling Nq t oc y/ Nm which is reproduced in fig. 
4d. The constant of proportionality depends on the initial 
occupation of the m = component. A larger initial oc- 
cupation number leads to a larger amount of atoms with 
a small kinetic energy. As a consequence more atoms re- 
main trapped reminiscent of the mechanism of evaporative 
cooling. 



In the present study the focus was set on a Id system 
but the results translate also to higher dimensions: For 
confinement in 2 and 3 dimensions ring [9, 12, 23] and 
shell-like [24] traps can be employed. Even more com- 
plex setups like arrays of BECs which provide a periodic 
'soft' or 'hard' trapping potential are conceivable. Such 
a scenario, which is reminiscent of the self-assembled lat- 
tices presented in Rcf. [25] can be realized via multi-well 
radio- frequency traps (see Ref. [24]). 
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